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Andreev reflection spectroscopy of ferromagnet/superconductor (FS) junctions is an important 
probe of spin polarization. We theoretically investigate spin-polarized transport in FS junctions in 
the presence of Rashba and Dresselhaus interfacial spin-orbit fields and show that Andreev reflection 
can be controlled by changing the magnetization orientation. We predict a giant in- and out-of-plane 
magnetoanisotropy of the junction conductance. If the ferromagnet is highly spin polarized—in the 
half-metal limit—the magnetoanisotropic Andreev reflection depends universally on the spin-orbit 
fields only. Our results show that Andreev reflection spectroscopy can be used for sensitive probing 
of interfacial spin-orbit fields in a FS junction. 


Spin-orbit coupling (SOC) is a key interaction in spin- 
tronics IDS!. allowing an electrical control of magnetiza¬ 
tion and, vice versa, a magnetic control of electrical cur¬ 
rent. In systems lacking space inversion symmetry- be 
it bulk, hybrid structures, junctions—SOC induces spin- 
orbit fields mil. as an emergent phenomenon. We are 
in particular concerned here with interfacial spin-orbit 
fields which are believed to be behind a wealth of new 
phenomena, not existent or fragile in the bulk, such as 
the tunneling anisotropic magnetoresistance (TAMR) [3J- 
[7 , interfacial spin-orbit torques [5] , or Skyrmions [9]. 

Interfacial spin-orbit fields are also important in 
semiconductor/superconductor [T0HI2] and ferromag¬ 
net/superconductor (FS) junctions [T3j for creating Ma- 
jorana quasiparticle states. It is the latter junctions 
that we focus on. We investigate the interplay of mag¬ 
netism and spin-orbit fields. We show that this inter¬ 
play leads to marked anisotropies in the junction con¬ 
ductance with respect to the orientation of magnetiza¬ 
tion. The most robust is the out-of-plane anisotropy 
(plane being the interface), which arises from the om¬ 
nipresent Rashba field [15]. A more subtle is the in¬ 
plane anisotropy, which arises from the interference be¬ 
tween the Rashba and Dresselhaus [TB fields, induced 
by a twofold anisotropy of the C 2 V type. A zinc-blende 
semiconductor (say, GaAs or InAs) as a barrier in an 
FS junction would create such an anisotropy, generat¬ 
ing spin-orbit fields C 2 V “butterfly” patterns, as shown 
by first-principles calculations m- Remarkably, the 
resulting magnetoconductance anisotropy—we term it 
magnetoanisotropic Andreev reflection (MAAR )—is gi¬ 
ant in comparison to TAMR, its normal-state counter¬ 
part, reaching a universal behavior in the half-metallic 
case. This is because Andreev reflection (AR) (which has 
no counterpart in the normal-state TAMR) is strongly 
influenced by interfacial spin-orbit fields. 

We specifically examine the influence of SOC and crys¬ 
talline anisotropy on the process of AR in which the re¬ 
flected particle carries the information about both the 
phase of the incident particle and the macroscopic phase 
of the superconductor to which a Cooper pair is being 
transferred m- AR is thus responsible for the proximity 


effect in which the phase correlations are introduced to 
a nonsuperconducting material [19l - [23] . While the main 
interest in AR is currently the proximity effect coupled 
with SOC, inducing Majorana states, in spintronics AR 
is used to probe the spin polarization in FS junctions Ba¬ 
sil. We argue that AR can also be a sensitive probe of 
interfacial spin-orbit fields. 

Our model FS junction consists of F (z < 0) and S 
(z > 0) semi-infinite regions separated by a flat interface 
at z = 0, with potential and SOC scattering. The scheme 
and possible scattering channels are illustrated in Fig. [I] 
For example, in conventional AR the incoming electron is 
reflected as a hole with the opposite spin, while spin-flip 
AR implies equal spin of the incoming and reflected par¬ 
ticles. These two AR processes, see Figs. 0 b ) and 0f), 
introduce, respectively, spin-singlet and spin-triplet su¬ 
perconducting correlations at the interface [22] 123] ■ 

We consider epitaxial-quality junctions, such as those 
used in TAMR [8], or point contact geometries [35], in 
which ballistic transport formalism is applicable. In 
diffusive tunnel junctions AR could be enhanced by 
electron-hole coherence [36]. In ferromagnetic junctions 
such effects would be absent for normal AR due to short 
coherence length, but spin-flip AR could be enhanced. 
(Ordinary effects of diffusion could be accounted for by 
renormalizing the tunneling parameters [32j ) . We gener¬ 
alize the Blonder-Tinkham-Klapwijk formalism m and 
solve the Bogoliubov-de Gennes equation (38, for quasi¬ 
particle states 'P(r) with energy E , 

(f; i ) 4,(r) = E * <r) - (i) 

with the single-particle Hamiltonian for electrons H e = 

-{H 2 /2)V Jl/m(z)] V- p,(z) - (A xc /2)Q(-z)m-a + H b ; 

for holes Hh = ~&yH*&y The unit magnetization vec¬ 
tor (see Fig. [l]) is m = [sin O cos <b, sin O sin <b, cos 0], a 
are Pauli matrices, A xc is the exchange spin splitting 
in the F region (Stoner model), m{z) is the effective 
mass, and /i(z) is the chemical potential. The interfacial 
scattering is modeled as Hb = (Vo d + w • <x) S(z), where 
V 0 and d are the barrier height and width, while w = 
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FIG. 1. (Color online) Top: FS junction. Magnetization vec¬ 
tor m is given by the polar angle 0 and azimuthal angle <f>. 
Current, /, flows perpendicular to the interface. To specify 
spin-orbit fields we use principal crystallographic orientations 
x = [110], y = [110], and a = [001]. Bottom: Scattering 
processes at the FS interface with SOC. Electrons (holes) are 
depicted by full (empty) circles. Vertical arrows denote the 
spin. The processes for a spin up incoming electron: (a) Spec¬ 
ular reflection, (b) Andreev reflection, (c) holelike transmis¬ 
sion, and (d) electronlike transmission, (e)-(h) Correspond¬ 
ing spin-flip counterparts. Spin-flip (equal electron and hole 
spins) Andreev reflection is in (f). 


[(a—/3)k y , — ( a+/3)k x , 0] is the effective SOC field includ¬ 
ing Rashba and Dresselhaus terms mm, parametrized 
by a and /?, respectively, for the crystallographic orien¬ 
tations see Fig. [T] The superconducting pair potential is 
given by A = AO(z)l 2 X 2 (the accuracy of such a step- 
function form of A is discussed in Ref. [39]), with the 
isotropic gap A. Similar methodology, for half-metal/S 
junctions with Rashba coupling inside the superconduc¬ 
tor was employed in Ref. [30]. With Rashba-only SOC 
one should still obtain out-of-plane magnetoanisotropy, 
and this is already implicit in this formalism [40] . 

Since the in-plane wave vector k|| is conserved, 
'F cr (r) = d' cr ( 2 )e lk l r H. The solution in the F region for 
incoming electrons with spin er is 
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with the spinors for the electronlike x% = (Xo-i0) T and 
holelike = (0,x_ CT ) quasiparticles, both containing 


X T a = ^oVl + (TcosOe l4> , vT^frcos©^ /V^2, (3) 

where a = 1(—1) corresponds to the spin parallel 
(antiparallel) to m. The electronlike (holelike) quasi¬ 
particle wave vectors in the F region are ka' h> = 
yJk‘ 2 F + 2 m F /h 2 [{-)E + <j A x . c /2] - 


In the S region the scattering states are 
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where the quasiparticle wave vectors are given by 
q e( - h ) = yjqp + (—)2 ms/h 2 \ZE 2 — A^ — fcjj. The super¬ 
conducting coherence factors satisfy u 2 = 1 — v 2 = 
(l + VE 2 - A 2 / E) /2. 

Using charge current conservation, the differential con¬ 
ductance at zero temperature, normalized by the Sharvin 
conductance [J Gsh = e 2 k 2 F A/ (27 rh) of a perfect contact, 
is 

G = E / [! + Ri(-eV) K(eV)] , (5) 

(7 F 


containing the probability amplitudes in the F region 


R^ h \E^) = Re(kl {h) 


e(h) 


2 -I- k e ( fe ) 2 
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which 


combine the coefficients for the scattering processes with 
and without spin flip for specular reflection and AR; V 
is the bias voltage and A is the interfacial area. 

To describe our results we introduce dimensionless 
quantities: Z = V^d^/mpms / ( h 2 ^/kpqp ) denotes the 
barrier strength EH EH, A q = 2 (Xyjm.pm.s/k? and 
\p = 2/3y/mpms/h 2 quantify the Rashba and Dressel¬ 
haus SOC, and P = (A xc /2) / jip defines the spin polar¬ 
ization in F. 

We first examine the influence of SOC on the FS con¬ 
ductance (see Fig.|2j) , for a metallic point contact (Z = 0) 
and for a moderate barrier (Z = 1). For the former case 
the conductance tends to decrease with increasing SOC. 
Even in the half-metallic case (P = 1) SOC does not 
give a finite subgap conductance; spin-flip AR is sup¬ 
pressed. In contrast, for moderate barrier (Z = 1), SOC 
enhances the conductance due to spin-flip AR, even for 
P = 1. Interestingly, at eV = A the conductance is not 
affected by SOC for any Z. Focusing on G(0), Fig. [2] 
shows that in a metallic contact increasing SOC steadily 
reduces G(0), while for a moderate barrier G(0) is a non¬ 
monotonic function of SOC, with a (P-dependent) max¬ 
imum which turns out to be due to spin-flip AR. 

The absence of spin-flip AR in metallic contacts can 
be explained analytically. For eV < A quasiparti¬ 
cle transmission is prohibited and subgap conductance 
G ~ Yla f ^ 2 kj| 2R^(— eV). In the half-metallic case 
the only contribution to AR comes from spin-flip AR, 
R„(E) = Re |r^ CT | 2 ^, because of the missing mi¬ 
nority spin subband in F. To lowest order in SOC and 
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FIG. 2. (Color online) Calculated normalized conductance, 
G(eV/ A), for different (indicated) spin polarizations P. 
Rashba SOC is \ a = 2 and Dresselhaus SOC is A^ = 0. 
Magnetization is in plane (0 = n/2). (a) No interfacial bar¬ 
rier (Z = 0), and (b) modest interfacial barrier (Z = 1) cases 
are shown. The dashed lines show G without SOC. The in¬ 
sets show the dependence of G( 0) on Rashba SOC. In (a), 
the subgap conductance for P — 1 vanishes for every \ a ,p', in 
the inset G(0) = 0 for this case. The illustrations summarize 
main qualitative impacts of SOC on conductance. 


Z then G(0) oc Z 2 Af with i £ {a,/3} [41] , vanishing if 
Z — 0. This perturbative quadratic dependence on the 
spin-orbit strength was also obtained in Ref. [10]. 

The calculated conductance features of SOC pH?[ - FR)| 
can be distinguished from k-independent spin-flip scat¬ 
tering by magnetic moments: For Z = 0 SOC always 
reduces the conductance and the subgap conductance 
vanishes for P = 1. In contrast, ^-independent spin- 
flip scattering m can increase the conductance and the 
subgap conductance is in general finite for P = 1. How¬ 
ever, similar features as those of SOC can arise in exotic 
superconductors without bulk inversion symmetry [48j . 

While the conductance changes are indicative of inter¬ 
facial SOC, magnetic anisotropy of the conductance is a 
true fingerprint. As the main contribution comes from 
AR, we call this anisotropy effect magnetoanisotropic 
Andreev reflection. We consider two configurations: in 
plane, in which magnetization m changes azimuthally 
(4>) in the interfacial plane, and out of plane, with polar 
(0) changes of m in a perpendicular plane (see Fig. 1). 
We define the in-plane MAAR as 


MAAR[ 110 ](d>) — 


G(0,O)-G(0,$) 

G(0,$) 


( 6 ) 


and the out-of-plane MAAR as 

• (7) 

$=- 90 ° 

The out-of-plane MAAR depends, in general on d>, but 
we choose the yz (<F = —90°) plane as its reference. 

The calculated MAAR, in Fig. [3j shows a nonmono¬ 
tonic dependence on SOC. For metallic contacts (Z = 0) 
MAAR is determined by the magnetoanisotropy of con¬ 
ventional AR. In the presence of a barrier (exemplified 
by Z = 1), MAAR gets strongly enhanced due to the ad¬ 
ditional contribution from spin-flip AR. In-plane MAAR 
exhibits C 21 , symmetry due to the interplay of Rashba 
and Dresselhaus fields, similarly to TAMR [2J d, 7\. If 
either of the two fields is absent, in-plane MAAR van¬ 
ishes. In contrast, out-of-plane MAAR is finite even with 
the Rashba field alone, which makes it a robust probe of 
this important interfacial SOC. Interestingly, at eV = A 
MAAR is always absent, as there are no effects of SOC 
on G here; see the discussion to Fig. 2. Additional effects 
(such as appearance of symmetry lobes) can arise due to 
the effective mass and Fermi wave vector mismatch [41] . 

Compared to TAMR , the magnitude of MAAR is gi¬ 
ant , varying by orders of magnitude upon changing the 
spin polarization P. (The experimentally measured in¬ 
plane TAMR in Fe/GaAs/Au junctions is less than a per- 


MAAR [1IO] (0) = 
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MAAR n 1Q] (0) [%] 


FIG. 3. (Color online) Top: Calculated in-plane magne¬ 
toanisotropic Andreev reflection (MAAR) with [110] crystal¬ 
lographic reference axis for Z = 0 (left) and Z = 1 (right) 
for different strengths of SOC A a = A^ = A at P = 0.4 and 
V = 0. Bottom: Out-of-plane MAAR with [110] crystallo¬ 
graphic reference axis. For Z = 0 the lines of A = 0.5 and 
A = 5.0 coincide. For the chosen reference axes and A a = A /3 
the in-plane and out-of-plane MAAR curves have the same 
magnitude and shape, but rotated to the corresponding ref¬ 
erence axis. 
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cent [8].) A detailed model comparison is shown in Fig. 
4 for both in- and out-of-plane configurations; TAMR is 
evaluated by setting A = 0. For a typical P of 40%, the 
ratio MAAR/TAMR is about 10. Moving towards half 
metals (P > 80%), this ratio climbs to more than 10 2 . 
This giant increase is best illustrated in the half-metallic 
limit of P = 1. For a weak SOC (which is typically the 
case) an analytical treatment gives [41] , 


MAARjho] ("£) 


2A Q A j a(l — cos 2$) 

Aa T A| + 2A a A/j cos 24* 


( 8 ) 


maar [ iTo] 


(A a + A/ 3 ) 2 (l — cos 20) 
3(A^TA|) —2A Q A / 3 + (A a + A j a) 2 cos 20 


(9) 


Therefore, the in-plane MAAR[ 110 ] (<f> = tt/2) ss 

4A Q A ( a/(A a — Xp) 2 , and out-of-plane MAAR[ lTo j(0 = 
7t/2) ss (Ac + Xp) 2 /(X a — Xp ) 2 , depending universally on 
the spin-orbit fields only , and diverging as A a ~ Xp (see 
the in-plane case in Fig. 4). In contrast, TAMR, which 
is proportional to the product X a Xp [7], has no singular 
behavior, and is not a universal function of A * only. 

We can trace this giant enhancement of MAAR over 
TAMR to spin-flip AR. Let us separate phenomenolog¬ 
ically the conductance G = G + G so into the sum of 
SOC independent and dependent parts. In TAMR typ¬ 
ically G^ G so , and TAMR ~ G so /G^ -C 1, even 



p 

FIG. 4. (Color online) Calculated in-plane (top) and out- 
of-plane (bottom) MAAR and TAMR as a function of spin 
polarization P for a moderate barrier (Z = 1) and V = 0. 
The in-plane case is calculated with X a = Xp = 1, while for 
out of plane we have included Rashba X a =0.5 only. The top 
inset shows the ratio of MAAR and TAMR for the in-plane 
case, while the bottom inset shows the color map of out-of¬ 
plane MAAR as a function of P and Rashba (or Dresselhaus) 
SOC Xi (where i could be either a or p). 


for P ss 1. But in FS junctions G^ 0 ) decreases with 
increasing P, eventually vanishing in the half-metallic 
limit. For P ss 1 the conductance of the FS junction 
is dominated by the spin-flip AR contribution to G so . 
Thus, SOC determines both the conductivity and the 
magnetoanisotropy. Furthermore, if A a « ±A,g, the spin- 
flip AR, and so the conductance, can be switched on and 
off by changing the orientation of m. For X a = Xp and 
$ = 0, m _L w and spin-flip AR yields a finite G. How¬ 
ever, if $ = 7t/2, then m || w and spin-flip processes 
are strongly suppressed; G(eV < A) at $ = tt/2 van¬ 
ishes. As a result, in-plane MAAR diverges if X a = Xp. 
Similarly for out-of-plane MAAR. 

There is one more peculiarity of MAAR in the half- 
metallic limit. If only Rashba (or only Dresselhaus) SOC 
is present, out-of-plane MAAR has a fixed universal mag¬ 
nitude of 100%. This is shown in Fig. 4 (in particular the 
inset for Xi < 1 shows MAAR of 100% for P « 1). It fol¬ 
lows from Eq. (|9]) that MAARj iT0 ] (0) ss (1 —cos 20)/(3+ 
cos20), which gives a universal amplitude of 100% for 
0 = 7t/2. In other words, G(0 = 0) = 2G(0 = tt/2). 
The origin of this universal behavior is traced to the spin- 
flip probability by scattering of spin-polarized electrons 
off spin-orbit Helds. The conductance is determined by 
spin-flip AR. For out-of-plane magnetization, 0 = 0, two 
fields, one along x and one along y, induce a spin flip. But 
for an in-plane magnetization, say along x, 0 = tt/2, only 
the spin-orbit field component along y can flip the spin. 
This gives the 2 : 1 ratio in conductances and 100% of 
MAAR. A more technical and detailed discussion of the 
differences between MAAR and TAMR can be found in 
Ref. SU 

Experimental realization of MAAR could follow the 
measurement geometry of TAMR Jj, ideally also the 
same junction, with the nonmagnetic metal that becomes 
superconducting at low temperatures. Magnetization of 
the ferromagnetic layer is typically rotated by an external 
magnetic field. This field can bring additional anisotropic 
orbital effects whose presence can be clearly identified 
from the field magnitude dependence |49j . However, one 
can avoid these extrinsic effects entirely if one uses dys¬ 
prosium magnets which can be oriented by the field, but 
do not need its presence to remain in the rotated posi¬ 
tion m■ Potential aspects of nonflat tunneling barriers 
can also be treated m ■ A practical alternative (espe¬ 
cially if ballistic junctions are desired) could be a point 
contact FS junction geometry (35, 52]. 

To conclude, we have applied a well-established theo¬ 
retical formalism to systematically explore the magnetic 
anisotropy of the conductance in FS junctions due to in¬ 
terfacial SOC. We predict a giant in- and out-of-plane 
MAAR—when compared with TAMR—exhibiting uni¬ 
versal characteristics in the half-metallic regime. The 
predicted magnetization control of the AR suggests a 
similar control of the superconducting proximity effect 
and Majorana states. Our findings reveal an unexplored 
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venue for AR spectroscopy, in the sensitive probing of 
interfacial SOC and related magnetoanisotropic phenom¬ 
ena. 
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TABLE SI. Dimensionless system parameters. 


Z 

barrier strength 

Vody/m F m s 
h 2 y/kpqp 

\ ol 

Rashba SOC 

‘loLy/mF'ms/h 2 

A /3 

Dresselhaus SOC 

2 fd^/mpms/h 2 

P 

spin polarization 

( A xc /2)/hf 

Fk 

Fermi wave vector mismatch 

qF/kp 

Fm 

effective mass mismatch 

m s /m F 


BOUNDARY CONDITIONS 

The reflection and transmission coefficients can be 
found by applying the following boundary conditions to 
the scattering states ensuring probability conservation 


v ]/ F \ = \p 5 | 

* a I2=0“ a I 2=0+ : 


(Si) 


JLA„n,s | = (w <t 

2m s dz 1 <t|2 = 0+ V 0 - 

+ {^n F l + V ° d ) 71 ^\^- 


(S2) 


where 


h 2 


(Aq \(f)ky 


W — — I — (Aq; + \p)k x 

2 y/m F m s \ 0 


(S3) 


represents the SOC field and // is 


V = 


1-2x2 


(S4) 


Equations (SI) and (S2) form a system of eight linear 
equations for determining the reflection and transmission 
coefficients. 


ANALYTICAL RESULTS FOR eV < A 

In the energy range below the gap (eV < A) transmis¬ 
sion of quasiparticles is prohibited. Therefore, the prob¬ 
ability coefficients for transmission are zero and proba¬ 
bility current conservation gives R% + R% = 1 so that the 
normalized subgap conductance can be written as 

O ~ V / dp dkkR^(-eV), (S5) 

^ a Jo Jo 

where we introduce a dimensionless wave vector magni¬ 
tude k = k^/kp; is the polar angle. The probability 
coefficient for AR is defined as 

R h AE, k) = Re (fcj \r h a J + k h _ a |r* i _ <r | a ), (S6) 

To derive analytical results we apply Andreev approx¬ 
imation, i. e., we neglect corrections of the order 
O (E / hf,s, A/ Hf,s)i with the chemical potential /j-f,s h 1 
the ferromagnet and in the superconductor, respectively. 
This leads to the wave vectors in F, k% ~ ~ kpk a , 

with the dimensionless wave vector 

k a = \/l + vP-k 2 . (S7) 

The superconducting wave vectors are q e ~ q h ~ kpq 
with the dimensionless wave vector 

q = fpl - *2. (S8) 

By inserting the scattering states to the boundary con¬ 
ditions [Eqs. (S1,S2)] we can obtain the scattering coef¬ 
ficients for AR with and without spin flip, respectively, 


/fc7 “ 16 Fhquv ( u 2 - v 2 ) s/F^ZW*_ a 
kp 


(S9) 


and 

r h 

' (7, — £7 


4 Fhquv 


( u 2 + v 2 ) q{k- a + k a ) + ( u 2 - v 2 ) [q 2 /F M + Fuk-.k, - 2 i y /F\ t F K Z(k- a - k a ) + 4 F K Z 2 } 


-4 Fhquv 


('u 2 - v 2 ) { [iVF M (k-a - k a ) - 4vdFVZ] W„, a + W% a + W a ,. a W^_ a } 


(S10) 
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with the denominator 

Q=(u 2 - v 2 ) 2 {q 1 + Fhklk 2 _ a + AF^F k Z 2 (k 2 + k 2 _„) + 16 F^F 2 Z 4 + 8 F K F M q 2 Z 2 } 

+ (u 4 - v 4 ) 2q (k a + k- a ) {F M q 2 + 4F K FhZ 2 + Fhk a k- a ) 

+ (u 2 + v 2 )~ Fl { q 2 (kl + k 2 _„) + ( u 4 + v 4 ) Fl,q 2 ik a k_ a 
+ W„,„ (u 4 - v 4 ) (k 2 _ a - k 2 a )UF% 2 q 

+ (u 2 - v 2 ) 2 (k-„ - K) { 2 iF^q 2 + 2 iF 7 J 2 k-„k a - 8 iF K F 5 J 2 Z 2 } 

+ (W 2 „ + W a ,- a Wl_ a ) {(u 4 - v 4 )F 2 M 2q(k a + k- a ) + (u 2 - v 2 f[-&F 2 M F K Z 2 + 2 F M q 2 + 2 F^k- a k„]} 
+ ( u 2 - v 2 ) 2 [wl„Fl (-k 2 _ a -kl + 2 k a k_ a ) + (Wl„ + W a ^W a ^ a W*_ a ) 2:«#’(*_„ - fc ff )} 


+ (w 2 ta + w^w; t _y 2 (u 2 - v 2 ) 2 f 2 m . (si 1 ) 

The dependence of the scattering coefBcients on the orientation of magnetization is included in the SOC matrix 
elements For in-plane magnetization (0 = n/2) they have the form 

W atCr = ak [(A q — \p) sin cos 4> — (A a + A^) cos yj sin <&] (S12) 

W a ,-„ = iak [(A Q — \p) sin (f sin $ + (A a + A^) cos y? cos 4>]. (S13) 

For out-of-plane magnetization and 4> = —7r/2 the SOC matrix elements read 

W a , a =j ak (Ac + Ap) sin ip sin 0 (S14) 

W a - a = k (A q + A^) sin y? cos 0 + iak (A a — A^) cosyp. (S15) 

The spin-flip AR coefficient in Eq. (S9) can be rewritten as 


r h a ^=aW;^[b+cW a , a +dWl a +eW a ^ a W^_ a +f{Wl a +W^W a ,. a W^)+g(Wl a +W a ^Wl_ a ) 2 ] @16) 


where all the coefficients are functions of k. 


ANALYTICAL RESULTS FOR HALF-METAL/S 
JUNCTION 

In the case of a half metal (P = 1) the only contri¬ 
bution to AR comes from spin-flip AR because of the 
missing minority spin band in the half metal. Therefore, 
the AR probability coefficient becomes 

/£(E) = Re[fc5|r*,| 2 ]. (S17) 

To derive analytical results for MAAR we neglect SPIN-ORBIT COUPLING EFFECTS ON THE 
higher order SOC terms. To lowest order in SOC [O (A?)] CONDUCTANCE 

the spin-flip AR coefficient [Eq. (S16)] reduces to, 

In the manuscript we show the conductance at a fixed 
r h a (^) }y* (S18) SOC for different spin polarizations P. Here we show the 

tT,<7 b(k) a ' a impact of SOC, in Fig. SI, on the conductance spectrum. 

For a metallic point contact (Z = 0) spin-flip AR is in- 
For weak SOC we can use Eqs. (S13,S17,S18) to ob- hibited. Increasing SOC decreases the conductance in F. 

tain the following approximate expression for the in-plane In half metals this leads to a zer0 subgap conductance, as 

MAAR for eV < A the on jy con t r ibution would come from spin-flip AR. For 

2A a A/3(l — cos 24)) finite Z and SOC, spin-flip AR is present and the conduc- 

MAAR[ 110 ] (4?) ~ _|_ A | _|_ 2X a Xg cos 24? tance exhibits a nonmonotonic dependence on the SOC 

strength. The conductance first increases with increasing 
Analogously, we derive an expression for the out-of-plane SOC and then decreases again (insets in Fig. SI). At the 


MAAR using Eq. (S15), 

MAAR - (A Q +A,) 2 (l-cos20) 

R[ii0](O) ~ 3( A 2 +A 2)_ 2AqA/S + ( Aa+A/3 )2 C0S 20 '* S 0) 

Remarkably, for small SOC the amplitude of MAAR 
in a half-metal/S junction is determined solely by the 
strength of the SOC parameters Aq, and \p. This is what 
we call in the manuscript the universal behavior. 
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FIG. SI. Calculated dependence of normalized conductance 
on bias voltage for different strengths of Rashba SOC X a and 
Dresselhaus SOC = 0 in an ideal system Fk = Fm = 1. 
The magnetization is in-plane (0 = 7 t/ 2). Panel (a) shows 
the results for a ferromagnetic system (P = 0.4) and panel 
(b) for a half-metallic ferromagnet (P = 1) in the absence 
of an interfacial barrier (Z = 0). In the insets a moderate 
interfacial barrier is present (Z = 1). 


gap edge, the conductance is not affected by SOC at any 
Z. 


MAAR VS TAMR 



FIG. S2. Left: Formation of side lobes due to finite mass 
mismatch, Fm = 5 (dashed dotted line) for Z = 1, P = 0.4, 
V = 0 and A Q = A /3 = 0.5. Solid line is for the no mismatch, 
Fm = 1 case. Right: MAAR as a function of the barrier 
strength Z. On the top is the in-plane MAAR (A a = A^ = 1), 
on the bottom is the out-of-plane MAAR (Aq. = l,A /3 = 0). 
In both cases spin-polarization is P = 0.4, V = 0, and there 
is no lattice mismatch. 


ent values of Z. This leads to a notable enhancement of 
the anisotropy in FS junctions with moderate Z with re¬ 
spect to TAMR. More remarkably, the MAAR amplitude 
can further be enhanced by orders of magnitude by in¬ 
creasing the spin polarization P, as discussed in detail in 
the manuscript. The experimentally measured in-plane 
TAMR in junctions with a single metallic F region is 
less than 1% [1]. In comparison to the measured TAMR 
values, the MAAR amplitudes at P > 50% can be con¬ 
sidered giant. 


The presence of superconductivity makes MAAR re¬ 
markably different from its normal-state analog of tun¬ 
neling anisotropic magnetoresistance, TAMR. The most 
important is the giant enhancement, as seen in Fig. 4 
of the manuscript. A more subtle difference appears in 
the symmetry pattern of MAAR. In the manuscript we 
show (see Fig. 3 of the manuscript) that the interplay of 
Rashba and Dresselhaus SOC induces characteristic C^v 
patterns of in-plane MAAR. If we also include a mis¬ 
match of the effective masses (say, if F is a ferromag¬ 
netic semiconductor), there can appear side lobes, while 
preserving the symmetry. An example for Fm = 5, in 
comparison with Fm = 1 is shown in Fig. S2(a). Side 
lobes may also appear in TAMR, but at much greater 
SOC magnitudes. 

A comparison between the barrier strength depen¬ 
dences of MAAR and the corresponding TAMR (A = 0) 
is shown in Fig. S2(b, c). The maximum value of MAAR 
and TAMR for a finite barrier occurs at slightly differ- 


MAAR VS TAMR: A PHENOMENOLOGICAL 
APPROACH 

As discussed above, MAAR and TAMR behave in qual¬ 
itatively different ways with respect to changes of system 
parameters such as SOC strength, barrier strength, or 
effective mass (Fermi wave vector) mismatch. To under¬ 
stand the physical origin of such differences we use a sim¬ 
ple phenomenological model proposed earlier for TAMR 
[2, 3]. Based on general symmetry considerations, this 
model predicts that, up to the second order in the SOC 
strength, the tunneling conductance has the form, 

G($)^G (0) +G' SO ($), (S21) 

where G^ is the SOC-independent part of the conduc¬ 
tance. The SOC-dependent part can be decomposed as 

G so ($) = Gj° + G“ iso ($), (S22) 
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where 


G‘ 8 °cx(A^ + A^) (S23) 

and 


G“ is °($) oc2A a A^cos(2$) (S24) 


are, respectively, the isotropic and anisotropic parts of 
G so . The in-plane TAMR is then given by, 


TAMR [110] ($) 


Gfg iso (o) - gr(^) 

G(°) + G so ($) 


(S25) 


In TAMR junctions interfacial SOC is a perturbative 
effect and the conductance is largely dominated by its 
SOC-independent part, i. e., G^ G so (<&). The TAMR 
can then be expanded as, 


TAMR [110| (<&) = (S26) 

Gams°(0)[i _ cos (2$)] G“ iso (0) cos(2$) 

G(°) f GW + "' ' 

Here we took into account that G™ lso ($) = 
G“ iso (0)cos(2$). Since G<°> > G“ iso , the in-plane 
TAMR is typically a small effect exhibiting an “8”- 
shaped angular dependence, 


TAMR[ 110 ] (4>) 


-cos(2$)] 

GW 


(S27) 


For the formation of side lobes higher order terms in the 
expansion in Eq. (S26) must be sizable and this occurs 
only for considerably large values of the SOC strength, 
for which GJq 1SO starts to be comparable to G^°). 

The situation is qualitatively different in the case of FS 
junctions, since here decreases when the spin polar¬ 
ization of the ferromagnetic lead increases. Thus, when 
the spin polarization is large enough and large effective 
mass mismatch [see Fig. S2(a)] is present, G^ can be 
comparable and even smaller than GJJJ lso . I n such a sit¬ 
uation higher order terms in the expansion in Eq. (S26) 
can become relevant even when SOC is weak. This ex¬ 
plains why in FS junctions with weak SOC the angular 
dependence of the MAAR exhibits side lobes, while they 
are absent in the TAMR. 

When the half-metallic limit is reached, the only mech¬ 
anism available for the subgap conductance is the spin- 
flip AR, which is determined by SOC. Therefore for 
half-metal/S junctions vanishes and the analog to 
Eq. (S25) reads as, 


MAAR[no] (4?) 


G an iso (Q) _ Ganiso W 

G so ($) 


(S28) 



FIG. S3. Color map of the magnitude of out-of-plane TAMR 
at 0 = 7r/2 for V = 0 and Z = 1 as a function of spin 
polarization and SOC strength A* with i E {a, (3}. 


Since the denominator in Eq. (S25) is determined by 
G(°\ the size of the in-plane TAMR is fundamentally 
limited by the value of the SOC-independent part of 
the tunneling conductance. On the contrary, the MAAR 
in half-met al/S junctions can in principle be largely en¬ 
hanced by making G so in the denominator of Eq. (S28) 
sufficiently small. As mentioned above, the only mecha¬ 
nism available for the conductance in half metals when 
the energy resides inside the superconducting gap is the 
spin-flip AR. Spin flip is however suppressed when the 
interfacial SOC field, w, and the magnetization orienta¬ 
tion, m, are collinear. Such a situation is not typical 
because the direction of the SOC is, in general, deter¬ 
mined by the carriers wave vectors [see Eq. (S3)]. Yet, 
when X a = Xp (A a = — A^) the SOC field aligns along the 
y || [110] (x || [110]) direction, independently of the di¬ 
rection of k [see Eq. (S3)]. Consequently, when A Q = Xp 
(X a = —Xp) and = 7r/2 (4> = 0), w || m, and the spin- 
flip AR is strongly suppressed. Under these conditions, 
the denominator in Eq. (S28) becomes vanishingly small 
and the MAAR[ 110 ](4> = 7r/2) is largely enhanced. This 
complements our discussion of the giant MAAR in the 
manuscript. 

In the out-of-plane configuration, magnetoanisotropy 
is also present if the system contains only Rashba or 
Dresselhaus SOC. As discussed in the manuscript, for 
moderate SOC (A* < 1) the out-of-plane MAAR in half 
metals is independent of the SOC strength and shows 
a universal amplitude of 100%. This universality is ab¬ 
sent for out-of-plane TAMR (Fig. S3). We can again use 
the phenomenological model to understand the reason 
for this. Analogously to Eqs. (S21-S24) we write the 
conductance as 


The approximate relation of the MAAR previously ob¬ 
tained in Eq. (S19) can now be reproduced by substitut¬ 
ing Eqs. (S23) and (S24) into (S28). 


G(9)«G (0) +G so (e) 


(S29) 
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FIG. S4. Color map of the ratio of out-of-plane 
MAAR/TAMR at 0 = 7r/2 for V = 0 and Z = 1 as a function 
of spin polarization and SOC strength A* with i £ {a:,/3}. 


and decompose the SOC-dependent part as 


G so (0) = Gj° + G“ iso (0), (S30) 

where 

Gif° oc 3A 2 (S31) 

and 

G“o iso (0) ex A 2 cos(20). (S32) 


A i with i £ {a:,/3} denotes Rashba or Dresselhaus SOC 
parameter, respectively. The out-of-plane TAMR is 


TAMR [1IO] (0) ! 


G“ iso (O)[l-cos(20)] 

GW 

Aj [1 ~ cos(20)] 

G(°) 


For the MAAR we get 

MAAR^yo] (©) = 


Gfg iB °(0) - G""(0) 
G so (0) 

1 — cos(20) 

3 + cos(20) 


(S33) 


(S34) 


The size of the out-of-plane TAMR is determined by the 
ratio of SOC strength Af and since G so (0)- 


Thus TAMR increases with increasing SOC. On the con¬ 
trary, the out-of-plane MAAR ratio in half metals is inde¬ 
pendent of SOC for weak SOC. Since only spin-flip AR 
contributes to the conductance, is suppressed and 
G so {©) becomes important. Therefore, both the denom¬ 
inator and numerator in Eq. (S34) are proportional to 
SOC parameters and the SOC dependence cancels out 
for a half-metal/S junction with either Rashba or Dres¬ 
selhaus SOC. This leads to a universal MAAR of 100% 
corresponding to G(@ = 0) = 2G(Q = 7t/ 2), i. e., for 
an incoming electron with out-of-plane spin the proba¬ 
bility for spin-flip is twice as large as for an electron with 
in-plane spin. If only one of both considered SOC param¬ 
eters is finite the SOC field lies in the x, y-plane and is 
isotropic with respect to <h. Hence, for out-of-plane spin 
in z direction two directions (x and y) exist, in which it 
can flip. Whereas, for an in-plane spin, e. g., in x direc¬ 
tion, only spin flip in one (y) direction is possible leading 
to the typical ratio of out-of-plane and in-plane spin-flip 
ratios of 2 : 1. Since in half metal the conductance is 
solely determined by spin-flip mechanism, this ratio can 
directly be verified in the MAAR amplitude. 


A fully numerical calculation of out-of-plane TAMR 
is presented in Fig. S3. Compared with Fig. 4. of 
the manuscript it is clear that MAAR is in general 
much greater. The corresponding ratio of out-of-plane 
MAAR/TAMR is in Fig. S4. The giant enhancement 
occurs for SOC parameters A^ < 1, that is, in the per¬ 
turbative regime, in accordance to our analytical calcu¬ 
lations. Particularly strong is the enhancement in the 
half-metallic regime, P > 0.8, in which out-of-plane 
MAAR is two orders of magnitude greater than out-of¬ 
plane TAMR. 
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